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SYNOPSIS 


G. MUKHOPADHYAY 
Ph.D. (Physics) 

Indian Institute of Technology Kanpur 
SEPTEMBER 1985 

INELASTIC ELECTRON SCATTERING FORM FACTORS 
FOR THE EXCITATION OF LOW-LYING LEVELS IN 
SOME 2p— If SHELL NUCLEI 


Inelastic electron scattering is a powerful method 
for investigating nuclear structure. By measuring the 
cross-section for the scattered electrons, one can determine 
the Fourier transforms of the nuclear charge and current 
densities between the initial and the final nuclear states. 

In comparison with ordinary gamma transition the momentum- 
transfer dependence of the nuclear matrix elements contains 
information about the spatial structure of the nuclear ground 
and excited states and can therefore give us a good test of 
the nuclear wave functions and the transition operators used 
in the theoretical models. 


Some recent inelastic electron scattering experiments 

4 - , - 

presented valuable data on the 0 2- as well as 

y 4 S * -L 

0 + -* it transitions in several 2r>- If shell nuclei, 

g.s. 2 , 

Within a measured momentum transfer range upto ^ 2 fm , 


the C2 form factors involving 0 


21 transitions are 
g.s. 1 ^ 


characterized by two peaks appearing at q Q ^£ 0.7 fm and 


-1 


at q ff ~ 1.7 fm" . Although the qualitative features of the 




XI 


form factors associated with the 0 + — — it transitions' are 

g.s. 2 

quite similar to those of the observed form factors for the 

+ 4. 0 

0 ► 1 transitions / the magnitudes of IF I for the 

Cj * s • X 

former transitions are smaller by an order of magnitude than 
those for the latter. 


in a recent shell-model study the C2 form factors for 
some Ti, cr and Fe isotopes were calculated in terms of 
wave functions resulting from empirical effective interactions 
operating in a restricted valence space involving the 

^7/2 + f 7/2 P 3/? con£ ^9 urat -^ ons • Tbe use of tile truncated 
shell model wavef unctions necessitated the use of highly 
mass-dependent effective charges for the interpretation of 
available data. 


In the present work we show that the angular momentum 
projected Hartree-Fock-Bogoliubov (PHFB) wave functions 
resulting from realistic effective interactions operating 
in the full 2p- If shell permit a consistent microscopic descri 
ption of the available experimental data on form factors in 
terms of nearly constant effective charges. To our knowledge 
the use of PHFB wave functions to interpret the available 
form factor data has not been reported so far. 

Chapter 1 presents the background and the genesis of 
the topics that we have selected for investigation in the 
present work. 

The choice of effective two body interaction is a 
crucial element of the • ca lcula tions involving HF or HFB wave 



XU 


functions, in Chapter n we have tested the various available 
effective interactions in the 2 p-lf shell in terms of their 
PHFB estimates for the subshell occupation numbers for 
protons and neutrons. The effective interaction providing 
the best agreement between the available and the calculated 
subshell occupation numbers is later used in the subsequent 
chapters for the detailed calculation of the inelastic 
electron scattering form factors. 


in Chapter III we present the calculation of the 


inelastic electron scattering form factors for the excitation 


of 2 ^ states in some even-even 2 p-lf shell nuclei in the 


PHFB framework using a slightly modified version of the kuo 
Brown effective interaction. Apart from the c2 form factors 
we have also calculated the static qnadrupole moments/ 


Q / and the reduced transition probabilities for e 2 transitions/ 


B(E2; 0 


+ 


> 2 ^), in these nuclei, in keen contrast with 


the results obtained in the shell model calculations 
_n 

1 7/2 

the use of the self-consistent PHFB wave functions leads not 


X 

involving f”/, + f 7 /? P 3/0 configurations, it is seen that 


only to a significant overall improvement in the calculated 
form factors, but also obviates the necessity of using widely 
different effective charges for various nuclei, in view of 
the intractibility of large-scale shell model calculation 
in the 2 p-lf shell nuclei* the method employed in the present 
work is expected to provide a reliable alternative for 
obtaining a reasonably accurate microscopic description of the 


xiii 


observed data on C2 form factors. 

in Chapter iv we present the calculation of the 

inelastic electron scattering form factors for the excitation 

of the 2 2 states in the even-even 2p-lf shell nuclei Ti 

and ^cr in the PHFB framework losing the modified Kuo-Bro\wn 

effective interaction. It is seen that the use of the yrast 

and yrare states with J = 0 , 2 projected from the self- 

consistent prolate and oblate intrinsic shapes permits a 

reasonably accurate description of the observed form factor 

data involving the 0 + * 2~ transitions. 

<3 .s . 2 

As mentioned earlier, an important feature that 

characterizes the observed form factors for the o + 

g.s. 2 

48 50 

transitions in the nuclei Ti and cr is their reduced 
magnitude * by roughly an order# compared to those for the 

+ *4 

O — 2-, transitions, present calculation offers a qualitative 
g.s. l 

4 - A. 

understanding of this feature of the ^ transitions 

in terms of the reduced overlaps between the initial (0 + ) 

and the final (2*^) states due to their different intrinsic 
parentage. 

Finally in chapter v we summarize the results and 
discuss some possible extensions of the present work. 



CHAPTER I 


INTRODUCTION 

The elastic as well as inelastic electron scatter- 
ing provides an excellent tool for exploring nuclear • 

1 2 

structure for a variety of reasons ' . The principal reason 
for the usefulness of the electron as a probe is that the 
interaction between the electron and the nucleus is the 
we 11- understood electromagnetic interaction that is descri- 
bed by an exact theory-quantum electrodynamics. As a 
consequence of this, the nuclear properties can be related 
to the electron scattering data in an unambiguous, quanti- 
tative way. Another feature that favours the choice of the 
electron as a probe is the weakness of the electromagnetic 
interaction; the coupling constant for the electromagnetic 
interaction is much smaller than the strength of the nuclear 
force. An electron, therefore, probes the nucleus gently, 
causing very little disturbance in the process. 

Considering the excitation of nuclear levels^ the 

| 

superiority of inelastic electron scat raring over other 
methods involving electromagnetic interact ions-and here 'we, 
may cite photoabsorption — lies in the separate variability 
of energy and momentum transfer. This allows us to map 
inelastic form factors for various nuclear excitations as 
functions of momentum transfer q. 



2 


Here we must also mention some of the inherent short- 
comings that are associated with the electron scattering 
experiments. Due to its small mass, the electron may be 
easily deflected/ and most of such deflections usually result 
in radiation. This leads to strong background to the exci- 
tation cross section. A theoretical drawback is also related 
to the fact that/ due to the strong distortion •£ the 
electron wavefunction by the Coulomb interaction, it becomes 
necessary to go beyond the plane-wave Born approximation 
(PWBA) for carrying out an analysis of inelastic scattering 
experiments involving heavy nuclei (2 > 20) . 

Within the framework of PWBA applied to the electro- 
excitation of low-lying nuclear levels of light nuclei, a 
measurement of the cross section for the scattered electron 
permits one to obtain the Fourier transforms of the nuclear 
charge and current densities between the initial and final 
states involved in the transition. The momentum-transfer 
dependence of the nuclear matrix elements contains infor- 
mation about the spatial structure of the nuclear ground 
and excited states and, therefore, provides us with a good 
test of the nuclear wavef unctions . 

A number of recent inelastic electron scattering 

experiments have provided valuable data on the 0 + -* 2 3 

3-7 • ■ 

transitions in several 2p-lf shell nuclei . Some of the 
salient features of the available form factor data covering 



3 


the range = 0, 1-2.0 fm 1 are the following; 

i) The form factors for the 0 + 2^ transitions are 

characterize by two peaks aopearing at q^^ 0.7, 

1.7 fm \ The form factors display a minimum at 
q ef £ - 1.3 fm" 1 . 

ii) The qualitative features of the form factors associated 

with the 0 + - 2* _ transitions are quite similar to 

+ + 

those of the observed form factors tor tne 0 -*■ 2^ 

transitions the magnitudes of tne form factors for 
tne former transitions are, however, smaller by an order 
of magnitude than those for the latter. 


The work presented in this thesis is motivated by a 

desire to provide a consistent microscopic description of 

3-7 

the available experimental data on the slectroexcitation 

of low-lying levels in the 2p-lf shell nuclei. This is 

sought to ba achieved by calculating the form factors in 

terms of the angular momentum projected Hartre e-Fock- 

8 

Bogoliubov (HF3) wave functions resulting from the realistic 
effective interactions operating in the valence space 

A 

calculation of the form factors in terms of exact shell 


spanned by the (lf 7 ^ 2 , 2p 3 ^ 2 , 2p 1/2 , orbits 


model wavef unctions generated within the ( If -y 2 , 2p 3 ^ 2 , 


n 


2pi/ 2 / lfj / 2 ) space is not tractable for the nuclei in the 
mass range A = 40-80 since the dimensionalities of the 
ralevant shell model matrices are prohibitively large. In 


4 


this context, the method outlined in the present work; is 
expected to provide a reliable alternative for obtaining a 
reasonably accurate microscopic description of the observed 
form factors. 

Variational wavef unctions for several valence 
particles occupying reasonably large configuration spaces 
usually depend quite sensitively on the effective inter- 
action employed in the calculation. We have, therefore, 
first examined, in chapter II, the suitability of various 
available effective interactions by comparing the ground 

state subshall occupation numbers resulting from these 

9 

interactions with the experiments . It turns out that a 
slightly modified version of the Kuo-Brown (KB) effective 
interaction^ 0 suggested by McGrory, Wildenthal, and Halbert^ 
provides reasonably good agreement between the available 
and the calculated subshell occupation numbers in several 
2p-lf shell nuclei. This effective interaction, hereafter 
labelled the MWH interaction, has been employed In subsequent 
Chapters for the calculation of inelastic electron- scattering 
form factors. In Chapter II we have also presented a compara- 
tive study of the projected HFB method as well as the spectral 

. 12 13 

distribution method developed by French and co-workars * 
vis-a-vis their predictions for the subshell occupation 
numbers. 

14 

We next discuss in Chapter III the calculation 
of the form factors for the 0 + -+ 2^ transitions in some doubly 



5 


even Ti, Cr and Fe isotopes in terms of projected HFB 

wavef unctions resulting f rom the MWH effective interaction. 

The quadrupole moments of the first excited 2 + states ag 

•well as the reduced transition probabilities for E2 transitions/ 

B(E2; 0 + - 2^)/ are also discussed in connection with the form 

7 

factors. Recently Iwamoto et a 1 . had carried out a shell 
model study of the 0 + -» 2^ form factors in terms of wave 

functions resulting from empirical effective interactions 
operating in the valence space of the ( If + ^ 7 / 2 ^? 3/2^ 
configurations. The calculation required proton and 
neutron effective charges possessing relatively strong 
number dependence for a reasonably successful interpretation 
of the available data. Our calculations involving an 
unrestricted valence space spanned by the ( lf^ 2 , 2p 
2p^jy 2 / lf^^ 2 ) n configurations demonstrate that the use of 
the self-consistent projected HFB wavef unctions obviates the 
necessity of employing widely different effective charges 
for various nuclei/ in addition to providing significant 
overall improvement in the calculated form factors* 


In Chapter IV we present the calculation of the 
inelastic electron scattering form factors associated with 
the 0 + - 2* transitions in the nuclei (^Ti, ^°Cr, ^Fe) 

J_ + 56 

and the 0 -*• 2^ transition in the nucleus Fe. The calcu- 

lations of the form factors involving the non-yrast 
levels employ the same overall calculatjjlonal framework as 
well as the set of input parameters - the effective fc^o-body 



6 


interaction as well as the spherical single-particle energies — . 

7 

as the one employed in Chapter III. Iwamoto et al . have 

4 - 4 - 

sometime ago also discussed the 0 -*• 2^ form factors in 

some doubly even Ti and Cr isotopes employing the configura- 
tion space ( f 7/2 + f 7/2 ^3/2^ con j unct i° n with the eirpirical 

1 6 

effective interactions given by Oda et al . as well as 

17 

by Yokoyama et al . . The earlier calculation brought out 
the sensitivity of the restricted shell model results for 
the 0 + - 2* form factors towards the choice of the effective 

interaction; the calculated results displayed large variation 
in going from one interaction to another. Also, none of 
the two interactions provided acceptable simultaneous agree- 
ment for the 0 + - 2+ as well as 0 + ** 2+ transitions. The 

calculations presented here demonstrate that a reasonably 

satisfactory interpretation of the available form factor 

n + 

data involving both the yrast as well as non-yrast J =2 
states can be provided in terms of the prolate and oblate 
variational states resulting from the MWH interaction. 


As pointed out earlier, the observed form factors 

for the o + - 2t transitions in the 2p-lf shell nuclei 
1/2 

are characterized by the property that [lF(q) l 2 (0 + -“2*)/ 

I F( g) | 2 (0 + - 2"tp] ~ 10. The calculations presented in 
Chapter IV provide a qualitative understanding of this 
feature of the 0 + 2* form factors in terms of the reduced 

overlap between the- initial and the final states due to 



7 


th ?ir differ-ant intrinsic parentage. 

Finally, in cnapter V we summarize the results 
obtained in, and wisdom gained from Cnapters II, III and 


IV 
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CHAPTER II 


PROTON— NEUTRON OCCUPANCIES OF SPHERICAL SHELLS IN 
THE GROUND STATES OF SOME.- DOUBLY EVEN NUCLEI 
WITH A = 46- VO AND THE EFFECTIVE INTERACTIONS 
OPERATIVE IN THE 2p-lf SHELL 

11*1 Introduction 

The expectation value of the subshell occupancy 
operator is one of the sinplest measurable dynamical quantities 
in the ground state of a nucleus. The measured values of the 
subshell occupation numbers provide a direct test of the 
nuclear wave functions involved in the microscopic description 
of the ground state* 

Recent years have witnessed considerable increase 

in our knowledge of the spectroscopic factors for pickup and 

i 1—6 

stripping reactions involving various 2p-lf shell nuclei 
The measured spectroscopic factors can be related to the sub- 
shell occupation numbers by using non-energy-weighted sum 
rules 7 . 

The calculations presented in this chapter are 
intended to examine the suitability of various effective 
interactions operative in the 2p-lf shell from the point of 
view of detailed spectroscopic calculations by comparing the 
ground state occupancies generated by these interactions vlth th< 
experimental values. 

Unlike in the case of the 2s- Id shell; a calcu- 
lation of the subshell occupation numbers in terms the exact 
shell model wave functions is not feasible in the context of 



the 2p-lf shell nuclei; the dimensionalities of the Hamiltonian 
matrices involved are prohibitively large. Furthermore, the 
collective nature of the observed properties of the nuclei with 
46 i a£ 70 necessitates the involvement of the full 2p-lf 
valence space, and any calculation errploying a truncation of 
the shell model space is not expected to give accurate results. 

The Hartree-Fock ( HF) or the Hartree-Fock-Bogoliubov 

( HFB) approaches are a priori quite suited for the calculation 

of subshell occupation numbers for the ground states resulting 

8 

from a given effective interaction . The requirement of the 
minimization of energy for the HF or HFB ansatz of the vari- 
ational wavefunction is connected primarily with the extraction 
of the one-body field. The lowest-energy HF or HFB intrinsic 
state, therefore, corresponds to the optimum single-particle 
field ( resulting from two-body interactions) appropriate for 
the ground states . This single-particle field in turn provides 
the best set of (single-particle) wavefunctions for calculating 
the expectation values of various one-body operators in the 
ground states. 

9 

Kota and Potbhare have sometime ago reported the 
calculation of the ground state energies and the subshell 
occupation numbers for several 2p-lf shell nuclei in the 
framework of the spectral distribution method developed by 
French and coworkers 10 /11 . The occupation numbers resulting 
from a nunber of effective interactions in the 2p-lf shell 


12 


■were compared. However, it turns out that the ground states 
in various 2p-lf shell nuclei are usually more than five 
widths (cr) away from the centroid of the energy distributions. 
This constrains severely the usefulness of the spectral 
distribution method involving just the first two moments of 
the Hamiltonian; the results obtained via the spectral distri- 
bution methods therefore become very sensitive towards the 
higher-order shape correction to the configuration level 
density. 

We have also presented in this chapter a comparison 
of the subshell occupation numbers resulting from the projected 
Hartree-Fock-Bogoliubov (PHFB) and the spectral distribution 
methods employing effective interactions with known spectros- 
copic properties; the comparison permits an assessment of 
the relative efficacy of the PHFB approach. 

Section II. 2 presents the calculational framework. 

A brief outline of the HFB method as well as the technique 
of angular momentum projection as applied to the calculation 
of the subshell occupation numbers are presented. In 
section II. 3 we discuss the theoretical results and their 
comparison with the available data. Finally section II. 4 
contains same concluding remarks. 
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IX . 2 Caleulational Framework 

11*2.1 The Hartree-Fock-Bogoliubov ( HFB) Method 

Consider the shell model Hamiltonian of the nucleus 
under consideration 

H =§<.alela> a J aa+(1/4)a 2 Y6 <aPlv A Iy6> a£a+a 6 a y ^ (II. 1) 

Where are the spherical single-particle energies and 
«xp 1 V A JY<5> is the antisymmetrized matrix element of an 
effective interaction. Here a a and a a are the creation and 
annihilation operators in the spherical basis la>. We wish 
to obtain a transformation from particle coordinates, to quasi- 
particle coordinates such that the quasiparticles are relati- 
vely weakly interacting: 


H 


E_ + H + H 


qp 


qp-int 


( II. 2) 


where E Q is the energy of the quasiparticle vacuum, H 
describes the elementary quasiparticle excitations, and 
H qp-int a ^ hopefully) weak interaction between the quasi- 
particles. 

The Hamiltonian (E + H ) may not, in general, 
preserve all the symmetries of H. in HFB one imposes cons- 
traints via the use of Lagrange multipliers such that some 
of the observables posse as desired expectation values. More 
specifically, we have 
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H' = H % - \ % , CH.3r 


where the Lagrange multipliers are chosen so that the number 
operators = 2 a^a^ and N y = 2 s.^ v a a v possess the expectation 
values 

< V % li o > = 2, <$ 0 !Nj,| § o > = A-Z (II. 4) 

In HFB theory one considers the general Bogoliubov 
transformation 





(u 


af3(3 


+ v ap a p ) 


( II. 5) 


Here U and V are NxN complex matrices in a basis spanned by 
N single-particle states. The 2Nx2N linear transformation 


t 
q i 


: u v "* 


4. 

a 1 ' 


= 

* * 





l v U J 


a 

V ^ 


is unitary. Writing 


M 



(II. 6) 


( II -7} 


the unitarity conditions MM ^ = M^M = I lead to the relations 


UU + W = UU + W = I, uv + vu = U V + V-U = 0, ( II. 8) 

+ <V, * 

where U , U and U denote, respectively, the adjoint, transpose 
and the complex conjugate of the operator U. 
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w* ^ *L 

Since M = M'ywe can invert the relations (II. 5): 


a 


(u pa q p 


+ w 


( II. 9) 


The quasiparticle vacuum is defined through the condition 


q a 1 > = 0 ( a11 a) 


(II. 10) 


Thus a solution to above equation is 


1 J A > = (normalization)* Ijq, v I0> 


•a 


( II. 11) 


where I 0 > is the particle vacuum*. 


Derivation of the HFB equations by the "equatlons-of-motlon" 

g 

method 

The density matrix P and the pairing tensor t are 
defined in terms of the expectation values of the operators 
ajijagand a^a^ , respectively , with respect to the quasiparticle 
vacuum given by equation (II. 10): 


V = <®o' a K ,S o >? t n t ) =< ®o ,a P a a r V • 

From the anticommutati«n relations for a^s, it follows that P 
is Hermitian and t is antisymmetric. One can easily evaluate P 
and t by using m""^ to convert a's into q' s y in conjunction 
with the condition (II. 10). Thus we have 
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p otp= < *o I a £ a a I *o> 

= <*o< (2 + qj + U 6a q 6 )l# 0 > 

= 2 l V Y P v 6 a 6 t 6* § V 6 P V - (vtv) ap - (II - 13 > 

( Here we have used the anticommutation relations 

C q a' q p> 6 aP'[ q a' q P = [ q a' ] = 0 ) • 

In a similar manner we can show that 

t = v+U * (11.14) 

As a next step, we apply Wick's theorem to replace the products 

of operators appearing in H' by the sum of normal products 

containing all possible contractions (A normal product 
*1* "t* "i* 

a.a _....*a * of particle operators is obtained by first writing 
<6 N * 

f 

these operators in terms of q 1 s and q‘ s, and then ordering 
them so that the creation operators are to the left of the 
annihilation operators. A sign (-1) P is included where p is 
the number of permutations to go from the original to normal 
ordering sequence) * 

We then obtain 

H* *= H' + H* + H* (11*15) 

O 2 4 i 

where involves n uncontracted operators! 
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H 


p = tr[(e-^ 71 N 7 t -\ y Ny + . |i> P+ |^it + ] f C ll« 16 ) 


H 2 ” ccfc (h %) a p * a + a a p ! 


+ lifi A ap :a H’ + ia S p ia a a p ! , 


C 11 . 17 ) 


H 


[ 4 = ( 1 / 4 ) 4 y6 <aPIV A lY6> ra+ a^ a 6 a Y s .. (II * 18) 


Here the HF Hamiltonian h, the HF potential p { and the pair 
potential is are given as 


h = e + r» Tap = Y g < ay I v A | P6 > , UX.19) 

% = l^ <a P |V A'V«> t Y 6 • < II - 20 > 

From the Hermiticity of p and the antisymmetric nature of t , 
one can immediately deduce that r is Hermitian and A is 
antisymmetric by invoking the properties of < IV A [ > { 

<apfv A lY6> = - <£a Fv a Iy6>; <qy f V A I p 6> = <p6lv A Ia Y > . 

The expectation values <$ I : al^a"^. . . ,a+ : I $> vanish 
by construction. Thus 

<§ 0 |Hi$ o >= tr [(e-V^/y^ ± T )P+ \ At+ ] . 

(XI. 21) 

Now assuming an "independent quasiparticle n form for the part 
we can write 

H 2 - i ! «’K 


( 11 . 22 ) 
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Equation (II* 22) leads to the commutator 

[ H 2 ' 4 1 = E a q a E a ? < + V «PV * (II * 23> 

On the other hand the use of the equation (11.17) results in 
the expressions 


[ V sS = s, [ ^yV^Y ^ V>P . 

(11.24) 


Equating the coefficients of a| and a^ we obtain the general 
HFB equations: 


/ h = h-^-X^UX.25) 


' h* a ' 


’5il 


r \ 

u i 

* i * 


-► 

= E 1 


v-A -h J 


V, 

t 1 j 


N 


where = (U^, U^,. . . .U^) . 


Time-reversal symmetry 

In our calculations we have divided all the basis 
states in the. configuration space spanned by the orbits 
A£ 7/2' 2p * 2p^/2 and lf 5 ^ 2 into two sets. The first set 
pontains the states I to, which are restricted to have 
(m + 1/2) equal to an even integer. The second set contains 
the time-reversed states ! k > = T!k>, which have (m + 1/2) = 
odd integer. The phase convention is T|nljm>= (-1) A ' + ^ m inl j-m>. 
The first set is 

^ lf 7/2,7/2* ' tl£ 7/2 / 3/2 ,lf 5/2 / 3/2' 2p 3/2 x 3/2 >/ ^ 1£ 7/2,-l/2*‘ 

^5/2,— 1/2' 2p 3/2,-l/2 7 2p l/2 / -l/2 > and Uf 7/2,-5/2' Af V2/-5/2^ 
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For quasiparticle operators 


< =p E< % a+ P + V ap a 5 > 


( 11.26) 




(11.27) 


we notice that p can not connect I k> and lk> . Further the 
only nonzero elements of t connect the states I k> to the 
states 1 k >. Since the interaction v conserves the magnetic 
projections, the potentials (r,A) are partitioned in the same 
manner as the densities (P, t) . 

Thus, 


P = 


p i 0 


0 P. 


/ t — 


0 t. 


S 0 


/ h = 


h x 0 


0 h. 


, A = 


0 A. 


A 2 ° 


(11.28) 


where P and h are Hermit ian. Since t and A are by definition 

CV* a <N> 

antisymmetric, t 2 = — t^ and A ^ = - A ^ . 

Substitution of (11.28) into (11.25) reveals that 
the energy matrix is also partitioned into two blocks and that 
the forms ( 11.26), ( 11.27) of the wavefunctions is retained. 

We have Imposed time-reversal symmetry by requiring that 


qT = T q+ T 


( 11.29) 


so that U Q - 

ap 


U aP' V ar ' V “P 


( XI. 30) 
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» 


The quasiparticle vacuum is now time— reversal invariant and 

P 2 = Pi r h 2 = h* , tj = t 1# A+ = A 1 . (11.31) 

Therefore only the block of the energy matrix that is related 
to (11.26) is diagonalized. 


Canonical representation for time-reversally symmetric 
HFB wavef unctions ( Bloch-Mess iah Theorem) 

We next consider the question of simultaneous diagonali- 
zation of p and t. We have just seen that these are Hermitian. 
Therefore we have to check their commutator. Consider 


A+ = 

r q*N 

/ 

B* = 

4 


q„ 

a 



a a 

v J 


Then A^" = M B^ • 


( 11 . 32 ) 


( 11 . 33 ) 


We can now define the generalized quasiparticle density 
matrix Q: 


Q aP =< S o' A p A a l *o > • (II - 34) 


Since <q^q 


a' 




<q* q|> = << 5, < 3p >= °' 


: a q p : 


, we have 


0 0) 

0 1 


( 11 . 35 ) 


The generalized particle density matrix is defined as 



< *, 



B a f § o > 


( 11 . 36 ) 
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Using the definitions (11.12) we get 


R 



t 

C'O 

1- P 


( 11.37) 


inverting relation (11.33) we find that R and Q are related by 
a unitary transformation: 

R = Mt Q M 


Since Q = Q, we have R" 


R. The latter relation yields 


P 



tt+ 


( 11.38) 


p t = t $ 


(11.39) 


Thus P^t^ = a nd this ensures that P-^ and t^ can be 

diagonalized simultaneously. One can therefore have a basis 

spanned by [ I k^> , Jk^> ; f k^> , 1 ic 2 >. . •••••] such that P 

and t have the canonical forms shown below: 



Here Jk> = • bj [ 0>and 



b 


t 

k 




b + = 2 C 
k a 


k/a 


a+ 

a 


(11.41) 



22 


Denoting the (real) eigenvalues of Hermitian and t^ by 
p k. 311(3 Hck res P ectiv ely , we obtain frctn the relation (11.37) 

1/2 

1 t _ 1 = [ P k (l“ p k ) ] • (11.42) 

2 2 2 

Setting = v k (u k + = 1) we get jt _j = U k V k # (11.43) 

kk 

Recalling the definitions (II. 12) we see that the quasiparticle 
vacuum can be expressed as 

I *o > - n + \ b k + b P I °> • ( 11.44) 


Approximations employed in the present calculations 

Inverting the special quasiparticle transformation 


t 

q k = 

* + ★ 

a k b k - v k b r 
k 

*A 

II 

! / b i + v k b k 

, (11.45) 

and substituting 

in ( II. 17) we 

get 



H 2 = 

H 11 + H 20 7 H 11 

■ & 

(K ll ) kk' “4 q k' 

• 

> 

H 20 = 

lik* ^ ( H 20^ kk‘ q t; 

+ 


. (11.46) 


As pointed out earlier/ H 2o acquires the form (11.22) when 

H* = 0 and Hi, is diagonal. Whereas the condition (H* 0 ) = O 

20 11 , 
leads to the BCS equations/ the requirements (H 2Q ) _ = O, 
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( H ll)kk' = 0 ^ Mk' ) lead to "HF-likS" equations provided 


h kk > = 0 and A 


kF 


O ( k / k' ) . We have made the latter 

approximation, viz. a = 0. The expansion coefficients 

kk T 

appearing in (11.41) have been obtained by diagonalizing 

in the spherical basis the HF-like potential h' which includes 

2 

the appropriate density P k = V k * 


h ap = < “ 1 e X V s v IP > + I <Slk IV a'P k > v k ■" . C I:E - 47) 

2 

The occupation probabilities V k are obtained by solving the 
BCS equations: 

(XI. 48) 




<kk ' V k ' k ' >l W 


The calculation involves iteration between equations 
(11.47) and (11.48) until a reasonable convergence is achieved 
in terms of both the expansion coefficients c | C/(X as we H as 


the occupation probabilities V k . Denoting by © k the eigen- 
values 
yields 


values of (11.47), it turns out that the condition ( H^) = 0 

/CrC ^ 


2e k u k v k ' A 


kK 


(u l - V?) = o 


( 11.49) 


The above equation leads to the following non— trivial implicit 

2 2 

solutions of (11.47) and (11.48) (here U k + = 1) : 


U 


1 n 

= 2 [ 1 + 






> < - - 


©, 




3 =; ] * (H.50) 
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Using (11,50) in (11,48) one gets 

A - = 3 S < k k [v A I k'k' > IsUSl . (11.51) 

A r~i — 

^ G k* + A k'k' 

2 

The condition 2 2 V, . = jsr vields 

£ k,Tc/ v 7i/y y 


5 t 1 

%/v 


©, 


]= N. 


/ 


_ 2 2 
S k + A kk 


71 A 


(11.52) 


II. 2. 2 Calculation of Subshell Occupation Numbers in Terms 
of States of Good Angular Momentum Projected from 
the HFB Intrinsic States 

Restricting ourselves to axially symmetric states I § > / 

we can label the states lk> by the expectation values of the 

operator j_* Thus the axially symmetric intrinsic deformed 
z 

HFB states with K = 0 given by equation (11.44) can be written 
as 


V = 


11 (uT + V™ bt bt ) I 0 > 
im 1 1 irn ii 


( 11.53) 


where the creation operators bt can be expressed as 


im 


irn 


S 

3 


(A. at ; b+ = 2 (_l) 

ii jm i 


im 


>m C Ji a+ j-n. - (II * 54> 


Here, the index i is used to 'distinguish between different 
states with the same m^ and j labels the spherical single 

particle orbitals lf 7/2 , 2p 3/2 , 2p 1/2 and lf^. The wave function 

7 12 
expressed in eqn.( 11.53) can be recast into the form 
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$o > = 


N exp 


a H> 1 0 > . UI - 


55 ) 


with 


ap 


2 

i 


C m « (ft 

j i V 

a r 


v*?a 

i iPi ■ i 


m 


U: 


a 


m a /_m p 


( 11 . 56 ) 


Here we have denoted the quantum numbers ( j a , by cc, 
and N is a normalisation constant. 

The states with good angular momenta J projected from 
the HFB state I § > can be written as 

Xx 


1 ^ > = 4 1 ^ K > =(^)^ K (Q) R(Q) l^dQ , (11.57) 

81? 

where R( Q) and cj^XQ) are the rotation operator and the rotation 
matrix respectively. 

In terms of the angular momentum projected states 
the expectation value of the operator fjj • 

~ + j r 


ri . = 2 a a , 

'j . jm im 

J m=-j J J 


( 11 . 58 ) 


can be expressed as 

nj = < r nj i 

= <*o' ( a jm V 4 1 *o> /<*o'4 1 *o > 

m=-j 

= / p(0) d J (©) Sin© 69/ 7 n(Q)6 J (9) Sin© d© 

o 00 Q ’ 

( 11 . 59 ) 


where 



The density matrix P(e) is defined as 


Here 


P (©) 


M(e) 


M(Q) 
1+M( ©) 


= F( ©) ft 


C 11.61) 

(11.62) 


and 


W s) 


J a J R 

S d .(©) d_ t (©) f . , . , 


(11.63) 


m'm 0 
a p 


Further we have 

n( ©) = v det ( 1 +m( ©) T (11.64) 

The subshell occupation numbers were calculated as 
follows. The results of the variational HFB calculation are 
summarized in terms of the amplitudes (U^ y V^) and the expansion 
coefficients These values are used to first set up the 

f matrix. We next evaluate F, M and l/( 1+M) for 20 values of 
the Gaussian quadrature points lying in the range (0,71/2). 
Finally , the subshell occupation numbers for the ground state 
are computed frcm equation (11.59). 

It should be mentioned here that an approximate 
expression of the subshell occupation numbers can easily be . 
obtained in terms of the expectation value of the operator T)j 
with respect to the intrinsic state I$ Q >, 

’Ij ~ < *o <*] '®o >= in lc ji' 2 (v ? )2 • lII - 65) 

13 

The expression (11.65) has in fact earlier been employed for 
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calculating subshell occupation numbers in seme doubly even 
Ni isotopes. It should however be emphasized that the estimates 
obtained from equation (n.65) provide the values averaged over 
the complete yrast band consisting of the states with J=0 to 
J = Jjnax . These 1 band- averaged ' values are lively 

to differ significantly from the estimates based on the use of 
the projected J=C wave functions given by Eq.( 11.59 )„ Tbis is 
because the considerations of angular momentum coupling require 
enhanced involvement of high-j subshells in the structure of 
high spin members of the yrast band. 


II. 3 Results and Discussion 


14 

In our calculations we have employed the Kuo-Brown 
(KB) effective interaction for the ( 2p 3/2' 2p l/2 / 1 ^5/2^ 

valence space as well as a slightly modified version of the 
KB interaction suggested by McGrory, Wildenthal and Halbert‘S 
(MWH). The KB interaction consists of the bare interaction 
derived from the empirical Hamada- Johns ton potential and the 
ranormal is ation due to the three— particle one— hole ( 3p— lh) 
excitation in the ^°Ca core (see Appendix A). McGrory et al . 
later found it necessary to modify the T = 1 part of the KB 
interaction with a view to optimize the agreement between the 
shell model and the experimental spectra for some calcium 
isotopes. The modified version of the KB interaction is 
hereafter labelled as the MWH interaction. The set of single- 
particle energies we have taken are ( in MeV) : e — 0.0/ 
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£( ' 2p 3/2 ) = 2 * 1 ' £(2p l/2 ) = 3 * 9 s* 16 e(lf 5/2 ) = 6 * 5 * 

11.3,1 Comparison of the observed subshell occupation numbers 
with the predictions resulting from the KB and the 
MWH interactions in the framework of the HFB method 
involving explicit angular momentum projection 

The results for the subshell occupation numbers for 
the lf 7 / 2 ' 2p 3 /2 ' 2p l /2 and lf 5/2 orbits obtained via the 
PHFB method for various Ti, Cr, Fa, Ni, Zn and Ge isotopes 
are given in Figure II. 1. The solid (dashed) straight lines 
join the theoretical estimates obtained with the MWH (KB) 
interaction. The experimental information about the spectros- 
copic factors for various pick-up and. stripping reactions is 
given in reference 1 for Ti isotopes , reference 2 for Cr 
isotopes, reference 3 for 56 Fe, reference 4 for Ni isotopes, 
reference 5 for Zn isotopes and reference 6 for ^°Ge. The 
experimental values for the sub-shell occupation numbers 

shewn here (filled circles) were extracted from the available 

g 

spectroscopic factors by Kota and Potbhare using the non- 

7 

energy-weighted French-MacFar lane sum -rules 1 . 

An examination of the results presented in Table 1 1,1 
and Figure II. 1 reveals the following; 

(1) The results obtained from the interactions KB 
and MWH'.: differ significantly in terms of the predictions 
concerning the (-r) orbit in the nuclei ' Cr, Fe 

2 % rV 

and 56/58 Ni. In the nucleus ^Ni, for example, the subshell 




FIG. II .1 COMPARISON OF THE OBSERVED SUBSHELL OCCUPATION NUMBERS 
WITH THE THEORETICAL ESTIMATES. 
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occupation numbers caning from KB and MWH interactions differ 

56 58 

by four particles. In the case of the nuclei * Ni even 
the predictions for the orbit l?2/2 by more than one 

particle. 


( 2) It is seen that botn the interactions KB and 
MWH predict identical results for the various orbits in cases 
of all the nuclei with A > 60. This convergence of the 
occupation numbers for the two interactions in the second half 
of the 2p— If shell is related to the fact that these inter- 
actions differ primarily in terms of the relative positions of 

7 7 

the centroids of the interactions operative in the a nd 

( ~ configurations. The relative importance of the orbits 
•L and -t, however/ decreases for the nuclei with A >_ 60 considered 
here. 


(3) In the nuclei 46 ' 4 ®Ti, ^°Cr and ^ Ni the 
predictions resulting from the MWH interaction are in distin- 
ctly better agreement with the experiments than those resulting 
from the KB interaction. In the nuclei 60 ' 62,64 Ni and °Ge 


the predictions of both the KB as well as the MVIH interactions 
are in reasonably good agreements with the experiments* In 
the nucleus 56 Fe, however/ the predictions resulting from the 
KB interaction for the various orbits separately, are in 
better agreement with the experiments compared to those 
resulting from the MWH interaction. This apparent superiority 
of the KB interaction is, however/ not sustained if one considers 
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the total strength for the 1 = 1 and 1=3 transfers. For 

example, the predictions of the KB interaction, viz, the 

values 2,90 and 7,13 for 1=1 and 1=3 respectively, differ 

significantly frcm the experimental values which are 1,93 and 

8.07 respectively. The MWH interaction on the other hand 

predicts values 1.54 and 8.46 which are within 20% of the 

experimental values. Interestingly enough, although it 

appears difficult at first sight to judge in favour of either 

of the two interactions in terms of their predictions for 

subshell occupation numbers for individual orbits in the 

54 

case of the nucleus Cr, a consideration of the total strength 
of the 1- transfers again brings out the superority of the 
MWK interaction. 


(4) A comparison of the theoretical occupancy 
estimates for the orbits 2p 3 / 2 ' 2p i/2 311(3 l£ 5/2 resultin 9 frctn 
the KB and the MWH interactions with the observed ones reveals 
a slight underestimation of the involvement of the f 5/ / 2 orbits * 
This may be implying the necessity of further modification of 
the part of the effective interaction MWH operating in the 
upper half of the 2P 3 /2' 2p l/2 an<3 ^5/2 s P ace * However, it 
is seen that the observed quantitative variation of the sub - * 


shell occupation number associated with the orbits 

2p and lfc/ 9 as a furction of the neutron number is repro- 

duced reasonably veil by the MWH- interaction* 
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11*3,2 Comparative study of the spectral distribution as well 
as the projected HFB methods 

In this section we compare the PHFB predictions with 
the estimates resulting from the spectral distribution method . 
The calculations are carried out with the KB as well as the 
MWH effective interactions. These calculations are intended 
to assess the relative efficacy of the two methods to bring 
out the differences in the two interactions in terms of their 
predictions for the occupation numbers. In Figure II* 2 we 
present our results for the orbits lf 7 y 2 3X1(3 2p 3/2* A corn P ar ^" 
son of the results obtained by the MWH- interaction with those 
resulting from the KB interaction brings out the relative 
insensitivity of the spectral distribution method towards the 
shell closure properties of the interaction. Compared to the 


KB interaction the interaction MWH' involves enhanced attra- 

2 

ction in the ( f 7 / 2 ^ states and an additional repulsion for 
the (£ 7/2 £ 3 / 2 ^ state s. Consequently/ the interaction MWH 
is expected to yield a near closure of the f^ 2 -shell for 
nuclei with N and/or Z = 28. The PHFB calculation successfully 
reflects this feature; the PHFB estimates for subshell occu- 
pation resulting from the two interactions in the nuclei ^ 2 Cr 
54 

and Fe display large differences. Although, the spectral 
distribution estimates are qualitatively consistent with the 
PHFB ones, one finds that the effects due to the subshell 
closure are significantly scaled down. It is not particularly 
surprising, since the involvement of a large number of single 




FIG. II .2 THE SUBSHELL OCCUPATION NUMBERS RESULTING FROM (a) THE SPECTRAL 
DISTRIBUTION METHOD AND (b) THE PHFB METHOD. 
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particle states a feature which excludes sub— shell 
closures — is an essential prerequisite for the applica- 
bility of the spectral distribution method involving just 
the first two moments of the eigenvalue distribution. 

^ II. 4 Conclusions 

A simple way to examine the suitability of various 
effective interactions from the point of view of detailed 
spectroscopic calculations is to compare the ground state 
occupation numbers generated from those interactions with 
the available experimental values. However, in view of the 
intractability of the exact shell-model calculations in the 
2p-lf shell it is worthwhile to employ simple yet reliable 
methods — such as the HFB method involving angular momentum 
projection — for the calculation of the ground state occu- 
pation numbers. 

From the results presented here it turns out that 
the MVTH interaction provides a consistently good overall 
description of the available subshell occupation data in a 
majority of the nuclei considered here. The results presented 
here also indicate that the spectral distribution method is 
not sensitive enough to permit a comparative study cf two 
effective interactions via only their occupancy estimates in 
cases where sub— shell closures are involved. 
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CHAPTER III 


INELASTIC ELECTRON SCATTERING FORN FACTORS FOR THE 
EXCITATION OF THE YRAST 2 + STATES 
IN SOME 2p-lf SHELL NUCLEI 


III.l Introduction 

In this Chapter we present the calculation of the 

Coulomb form factors for the 0 + - 2^ transitions in the 

doubly even nuclei 46/ ^ r8/ ^°Ti / ^°'^Cr and ^'^Fe. Apart 

from the Coulomb form factors we shall also discuss here 

the results of the calculations’ of the static quadrupole 

moments , Q , for the 2"t" states, and the reduced transition 
2 + 1 

probabilities for E2 transitions, B(E2; 0 + -* 2 ^) , in these 
nuclei. 


Inelastic electron scattering provides a powerful 

1 2 

method for investigating nuclear structure. ' This is 
Decause the momentum— transfer dependence of the nuclear 
matrix elements contains information about the spatial 
structure of the ground and the excited states, in conpari- 
son with the ordinary gamrna transitions. The calculations 
presented here demonstrate how the available form factor 
data can serve to characterize the essential structural 
differences in various nuclei in the A = 40-60 region. 

Recent years have witnessed a significant, increase 
in the experimental data on the 0 + - 2+ transitions in the 

2p_lf shell nuclei 3 * 7 , within a treasured momentum-transfer 
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range upto <3 e j=f~ 2.0 £m , the C2 form factors are character-?,, 
ized by two peaks appearing at q Qf f ^ 0.7 fm’ 1 as well as at 
^ * The magnitudes of these peaks are expected 

to provide a sensitive test of nuclear models* 


n 

In a recent she 11- model study by- Iwamoto et al . , the 
C2 form factors for some Ti, Cr and Fe isotopes were calcu- 
lated in terms of wave functions resulting from empirical 


effective interactions operating in restricted valence soace 
involving the f^ 2 + f^” 2 P 3 / 2 configurations. Their 
calculations , which involved proton and neutron effective 


charges possessing relatively strong number dependence, 

showed that the systematic deficiencies of the f^ 2 model, 

vis-a-vis the relative magnitudes of form factors at the 

n 1 

two peaks, can be partially remedied by the + i"/2 P 3/2 

model. In this cnapter we show that the projected Hart re e- 
Fock-3ogoliubov (PHF3) wave functions resulting from appropriate 
realistic effective interactions operating in full 2p-lf 
snell permit a consistent microscopic description of the 
available data on C2 form factors in terms of nearly constant 
effective charges. 


In Section III. 2 we discuss some details of the 
calculat ional frame work. Section III, 3 discusses the 
comparison of the calculated Coulomb form factors and E2 
transition strengths involving the 0 + and 2^ levels, as 
well as the static quadrupola moments for the 2* states 
with the available experimental results. Section III. 4 
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contains some concluding remarks. 


i;fl.2 Calculat iona 1 Framework 

III. 2.1 Inelastic Scattering to Discrete Nuclear Levels 

2 

In this section we present an outline of the 
basic formalism relsvant for a discussion of the electro- 
excitation of nuclear levels. 

We consider tne inelastic electron scattering in 

the one-photon-exchange approximation. Let an incident 

electron with four-momentum k^ = (k^, iE^) be scattered 

— ► 

through an angle Q to the four-momentum k 2 ^ = ^2' ^2^’ 

The process involves the exchange vjf a Virtual photon with 
four— momentum u = (q, iE) . For electron energies of 
interest (hundreds of MeV's) wa can take E^ ~ k^ and = ^ 2 * 
The process of Inelastic scattering satisfies the kinematic 
equations 


ki-k 2 


= q 


= ( k, + k" 


- 2k- L k 2 cos©) 


1/2 


(III.l) 


k 1 + M = k 2 + E' 


with 

B' 



, 2 . 1/2 

+ M ) 7 


(III. 2) 


( III. 3) 


Here M and M' are the target masses before and after exci- 
tation, respectively. 

The Hamiltonian describing the interaction of the 
nucleus with the electromagnetic field (X, i$ e ) of the 
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electron is 

H = H c + r^r , (III. 4) 

where the Coulomb part is 

H C = ® 1 d3r (III. 5) 

and the transverse part is 

= -e S [ j c ( r) # A(r) + jt s ( r) .vx A(r)] d 3 r „ ( III. 6) 

Here eP^(r) Is the nuclear charge density ( with SP ( r)d 3 r=Z) y 

e j (jt) is the nuclear current density and ejl ( r) is the 
c s 

nuclear magnetization density, (The actual form of appropriate 
quantum mechanical operators relevant for the inelastic 
scattering will be given a little later in this section). 

From now onwards we consider only the Coulomb Inter- 
action. Writing § ( r) as 

e 

§ (r) = /[P e (r')/fr - r' I ]d 3 r' , (III. 7) 

we obtain 

H c = e 7[P n ( r) P Q ( r' )/| r - r'|]d 3 r d 3 r* # (III. 8) 

— 1 -♦ ^ 

Expressing R (R = |r - r'l ) as 

( 1/R) = ( 2 "Jt ) ” ^/exp ( iq.&) [e dg * R /R'] d 3 R' d 3 q 

and evaluating the integral w.r.t. R' we get 
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Putting in plane waves to lowest order for the incident and 
outgoing electrons, 

y ^ = u( Joj.) exp i(k 1 *r-ot) / '? f = u( i? 2 > exp i(ic 2 *r -cot) / 

( 111*13) 

and enploying the expression = for the transition 

charge density, we can simplify the expressioi.3 (III.-ll)j 

-+ -► -► -♦ 

X(g) =«e / u( k 2 ) exp iC-k^.r'+cot) Y 4 expi( k^r'-wt) uOc^ 

t ^ 3 

exp ( - i q.r 1 )d r* 

= -e ( 2it) 3 6( k^ic^q) (u(k 2 ) Y 4 uCi^)) . (III. 14) 
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From equation (III. 10) we then obtain 


H C = l6t 2 e / q 2) (u(k; 2 ) y ^u( k ] ,) ) 


S 

LM 


Y LM (q) f m (g) ■ 


( III. 15) 

In the framework of tne one-pnoton exchange approxi- 
mation, the differential cross section is next obtained by 

2 

multiplying IH I by 2ft times the appropriate phase space. 

The latter is given by 


( 2ft) 3 d 3 k 2 6(1^ + M - k 2 - E') ; (2=1) . (III. 16) 

2 

As pointed out by Uberall , the dynamic (nuclear) recoil 
effects, involving the free motion of the nuclear centre of 
mass, do not affect at all this phase-space factor which 
summarizes just the kinematic recoil effects. The cross 
section in the laboratory frame is thus 

d a = ( 27® 1 A 1 ) 6(k 1+ M-k 2 -B')(m 3 /E 1 S 2 ) 1 H c l 2 d 3 k 2 /( 2*) 3 . 

(III. 17) 

Here the factor (m 3 /E 1 E 2 ) is due to the invariant normalization 
of the spinors u(k 2 2 ) and ( k-j/E^J is just the incident electron 
f lux. 

Considering the scattering of the outgoing electron 
into the solid angle dQ we obtain 
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S 6(k 1+ M-k 2 --E'(k 2 )) d 3 k 2 = k 2 dQ ll+dE'(tc 2 )/dk 2 r 1 

= k 2 dQ[l+(2k 1 /M)sin 2 C©/2)]" 1 / 

(III. 18) 

since we have dE'(k 2 )/dk 2 = (dE’/dq) (dq/dk 2 ) 

= (q/E 1 ) (dq/dk 2 ), (III. 19) 

using the equations (III. 3) and (III.l) and employing the 
approximations E' = M, k 2 = k^. 

2 

The term I HI, I in equation (III. 17) involves the 
expression lu(k 2 ) Y 4 u( k^) 1 . The actual cross section 
involves a sum of this expression over final spin states and 
'an average over initial states; it should therefore be 
replaced by 


1 

2 


s/, 3l ru <*2' S f )Y 4 u(k l' S i>' 


The spin sum can be rewritten as 


( Hi* 20) 


= Sl V K 2' s f> (VafsVVh 1 “6 lE l / s l >CY 4 , 4cr u <; e 2' S f ) • 

( III. 21) 


It can be snown that 

I 1 u p ( h' S i> “5 ( ''l' S i ) = [(X 1 + ira e ) / 21I "Jp6- 


Tnus the expression (III. 21) simplifies to 
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Tr [( (X 2 +im e )/2im e ) Y 4 ( (X 1 +im e ) /2im e )y ^ j 
= ( E i E 2 + + m e' ) / m e * (III. 22) 


We next replace the nuclear matrix element f_..(q) by 
its operator definition/ viz. 


f LM ( 5) = i L / P op ( j L (qr) d 3 r 

= iL i e k J L (9 r k ) ” here P ° P (i>= l e k 6 (i “ ? k ) ' 

( III. 23) 

and consider the transition matrix elements-CO^MJ (q) 1 J.M.>/ 

f f IM ^ xi 

where J^M^(J^M^) denote the initial (final) nuclear states. 

We then obtain the differential cross section, expression 
(III. 17), by using the equations (III. 15), (III. 18), (III. 22) 
and (III. 23). Starting out -with the ground state 
the operators hJ( H^) in I H^J can lead to the final states 
lj £ M f > ^ £ M£>). The differential cross section can, there- 

fore, be looked upon as a "density matrix" with indices (J^M^), 
(J£M£) : 


(dq/dfi) 


J i H> J f ' J f 
M i^ M f' M f 


= (471) 2 (2 
2 


( J i M i ,LMlj f M f )( J^, 


X R. 


.MM’ 


LL' 


L'M* I J f M f ) 
( 111^24) 


wnere. 


MM* 

^LL' 


= v e> 


g*' 


if* if 
Y IM (q) (f L > f L- 


(III. 25) 
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v c (e) 


^ 2 /q 2 ) 2 (2^2 + k 1 .k 2 + rtig) 


( III. 26) 


and 


if, 


£ L^ (( 3 ) =< JfH f° P (q)l ! = [J f ]<J f M f I f°*J 1 J ± *j>A 

(Here [a] = ( 2a + 1) 1/2 ) # 


( III. 27) 


Finally t by setting = j , = M f and summing over 

final spin projections and averaging over the initial ones 
we obtain the differential cross section for electro- 


excitation ( Jj_-+ J f ) i 


d a 


J i - J f 


dQ 


) 


Coulomb 


r j- 2 „ ,d0v J i ** J f' J f 

“ [J J M^M f Mf/Mf - (III - 28) 


The sums over the spin projections can be evaluated 
by invoking the orthogonality of the Clebsch— Gordan coefficients 


2 -1 -1 


M 


^ M ^(J i M i ,IJ4|J f M f )(J i M i/ L , M , l J f M f ) =[J f ] [L] [L ( ] 


6 ll , 6 mm‘ ' 


and the formula 


1 * m (q> = MV** . 


The expression (III. 24) thus yields 


J j -*■ J ^ 

d<? / ^2_2, *4 % v .2 


). ( 22 z e^ 2 /k 1 A'*) v c (0) I F(q) I £ r0ooil ,( HI. 29) 


where 


' reco 11 


= [1 + ( 2k 1 /M) sin 2 (©/2)] 


-1 


( m*30) 
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2 

and the form factor lF(q)I is defined as 

I i'( q) 1 2 = (4VZ 2 ) ( 2J.+1)" 1 S Kj II f° P (q)lj >l 2 . 

X J_j £ Xj X 

(III. 31) 


III. 2. 2 Inelastic Scattering Form Factors in Terms of 

Projected Hartree-Fock-Bogoliubov wavef unctions 

The framework of angular momentum projection 
has already been given in some detail in the preceding 
'Cnapter. Here we shall just outline the method as it applies 
to the calculation of the C2 form factors. It may be pointed 
out that, although the PHFB metnod has earlier been success- 
fully employed for the calculation of the level energies, static 
quadrupole moments and reduced transition probabilities for 
E2 transitions, its use in the calculation of C2 form factors 
has not been reported so far. 

Using tne projected HFB wave functions [see 

equation (11.57)], one obtains 


< TP 


0 


! f 


op 

LM=0 



> = 






x 


* <3.^(9) n(9) { ^e^nlj^qr^lnp 1 p > x 
x<alYyPip> J slnede . (III. 32) 

The density matrix P( 9) is defined by equation (11.61).. 
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Here the normalisations 



are given as 


J 


n 


V2 

X 

0 


n( ©) 



(9) sin9d© 


( III. 33) 


where 

n( 9) = ^det (1+M(e))' 


( III. 34) 


The reduced transition probability for electric 

quadrupole transition, B(e 2, 0 + -* 2t) , and the static 

quadrupole moments for the first 2t state, Q , are given 

1 2 + 
by 


B(E2, 0 + - 2+) = 


iH ' < *0 11 q2|is o> |2 


( III. 35) 


and 


' 2 + 


<*2 ,Q 0 > * 2 > 


(III. 36) 


where, 

= 7^ (^ 5 ) (Q) (III. 37) 

and b is the oscillator length parameter. 

The expression for the reduced matrix element appearing in 
eq.( 111,35) can be obtained from eq.(III.32) by replacing 
j L (qr) by r L . 

Tne calculations were performed as follows. First, 
the wave functions l$ 0 > were obtained by carrying out the 
HFB (or the HF-BCS) calculation. The details of the intrinsic 
states-the amplitudes (U™, vj) and the expansion coefficients 
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m , 

C - are used to set up the f matrix (see aquation (11.56)). 

ji 

Then F, M and l/( 1+M) appearing in aquations ( II. 61) -( 11.63) , 
are evaluated for twenty values of the Gaussian quadrature 
points between the range (0,it/2). Finally, the form factor 
is computed from equation (III. 32). 


For the evaluation of the sing la— particle matrix 

elements, we have employed the oscillator wave functions with 

the length parameter given by b = 1.01 x A 1 / 6 fm. The center- 

1 2 

of-mass correction ' has been taken into account by multi- 

2 2 2 

plying a factor exp(b q /4A) . The correction due to the 

finite size of proton has been incorporated through a factor 
2 2 

axp(-a q / 4A) where a = 0.59 fm. 

ir Jr 


Further, in order to compare the form factor calct*- 


lated by PWBA with the experiment, the experimental data 

7 

have been plotted at the effective momentum transfer <3 e ££ 
instead of the kinematic q: 

2 


*eff 


r , 3VT Ze" 

- ^ [i + 27? W ] 


( III. 38) 


where R is the rms radius and E is the energy of the incident 
e lectrons. 


III. 2. 3 The Input Parameters of the calculation 

In our calculations we nave enployed tne MWH effective 

Q If) 

interaction * described in the preceding Chapter. As 
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demonstrated earlier/ this effective interaction is quite 

appropriate from the point of view of the subsnell occupa— 

tion numbers as well as the deformation trends 11 in the 

2p-lf shell. The single-particle energies for tne 2p-lf 

41 

shell orbits are taken from the Ca spectrum: 

e (lf 7 y 2 ) = 0.0, e(2p 3 ^ 2 ) = 2.1 MeV, £(2p = 3.9 MeV, 
and £(1£ 5 ^ 2 ) = 6.5 MeV. 

III. 2.4 The Intrinsic States for N=28 and N^28 Isotopes 
We have obtained the axially symmetric intrinsic 
states for the N^28 nuclei 46/4 ^Ti, 5 °Cr and ^ 6 Fe by following 
the HFB procedure discussed in Refs. 12-14. Pairing corre- 
lations between only the like particles are allowed in our 

50 52 54 

calculations. For the N=28 nuclei Ti, Cr and Fe, however, 
we have employed the deformed HF-BCS procedure discussed by 
Goodman 1 ^. It turns out that, whereas the HFB prescription 
leads to intrinsic states possessing sizable quadrupole 
deformations in the case of the N^28 nuclei, it yields nearly 
spherical intrinsic states in the nuclei with W=28 closed shell. 
This is because the HFB framework, which is expected to treat 
the deformation and the pairing effects on the same footing, 
overemphasizes the latter whenever a subshell closure is 
possible. The near— sphericity of the N=23 HFB states implies 
vanishingly small (<10%) amplitudes for angular momentum 
states with jyb, a feature which renders them unsuitable in the 
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context of the present calculation which requires a satis- 
factory description of both the J=0 as well as J=2 states. 
The HF-BCS intrinsic states for the i\f=28 nuclei, on the 
other hand, are quite deformed. 


III. 3 Results and Discussion 


Since the J = o' 1 ", 2 + states projected from the HFB 

intrinsic wave function consist dominantly of (lf 7 , 2 , 2p 3 ^ 2 ) n 
configurations, the matrix elements <flj 2 (qr)[f> and 
<fl j 2 (qr) I p > are expected to play important roles vis-a-vis 
the observed squared form factors; in fact the gross quali- 
tative features, viz. the appearance of two maxima at q N 0.7 
fm ^ and q <s * 1,7 fm \ can be understood in terms of the 
q-dependence of these matrix elements. With the assumption 
that the q-dependence of these matrix elements is not altered 
significantly by the effective charge, it turns out that 

(see Fig. Ill, 1), whereas the matrix element I <£ I j 2 ( qr) I f > j 

-1 -l 

possesses maxima at q ~0.8 fm and q ^ 1.9 fm , the maxima 

of the matrix element I <£ lj 2 (qr) | p>J occur at q <^0.7 fm -1 

and q ~ 1.6 fm" . Tne zeros of the matrix elements <f Jj 2 (qr)lf> 

and <f I j 2 ( qr) | p> occur at q ~1. 6- fm ^ and q o»l. 2 fm 1 , 

respectively. Since F(q) involves (see Eq.(III.32)) a linear 

combination *f the matrix elements <flj 2 (qr) I f > and 

» 

<£ | j ( qr) |p >, it is expected to exhibit two maxima in the 
ranges q = (0.7 , 0.8) fm and q = (1*6, 1.9) fm , as well as 
a zero in the range q = (1.2, 1.6) fm . 


This is consistent 




FtG.IIl .1 THE MOMENTUM-TRANSFER DEPENDENCE OF THE MATRIX 
ELEMENTS OF jjlqr) IN THE OSCILLATOR BASIS- 
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2 

with the observed features of lF(q) 1 in various 2p-lf 
snell nuclei* 

III. 3.1 N=2b Isotones 

2 + + 

The squared form factors [ F( q) | for the 0 -» 2^ 

transitions in the N=28 isotones are calculated and compared 

with the experiments in Fig. III. 2. 

In our calculations we employed the effective charges 
= 1.7e for protons and e^ = 0.7e. for neutrons for the 
three nuclei 5 °Ti, ^ 2 Cr and 54 Fe. As we shall see later, this 
set of effective charges is consistent with the available 
electromagnetic properties involving the 0 and 2^ states 
of these nuclei. From the results presented in Fig. III. 2, 
one finds that the overall quantitative agreement of the 
calculated form factors with the experiments is reasonably 
good, particularly in view of the fact that the effective 
charges were not fine-tuned for each isotone. 

A noticeable discrepancy, however, occurs in the 
case of the nucleus 54 Fe for q e££ > 1 * 4 fm ' the present 
calculation underestimates the magnitude of the form factor 
at the second peak; around q e££ ~ 1 * 7 ±m bir about 34 percent. 
This may be due to the inadequacy of PWBA for large q Qff and/or,' 
insufficient configurational admixtures in the HF-3CS wave 
functions in 54 Fe; the latter feature is also reflected in 
the low value of the contribution of neutrons towards the total 
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54 

intrinsic guadrupole moment for Fe (see column 2, Table 
III.l) . 


We have also given in Fig. III. 2 the results obtained 

7 

by Twamoto e t a 1 . in the framework of shell model involving 

the 2 + f 7/2 P 3/2 GOnf i 5 ur a ti°ns. The effective charges 

(e , e ) = (1.4e, l.Oe) , (l„6e, 1.2e) airi (1.8e, 1.4e) were 

50 52 54 

employed for the nuclei Ti, Cr and Fe, respectively, 
in these calculations. 


It is seen that .whereas the PHFB results are quite 

5 4 

close to the shell model ones in the nucleus Fe, they show 

distinct improvement over the latter in case of the nucleus 

52 Cr. The shell model calculations, however, yield better 

agreement with the experiments than tne PHFB ones in case of 

the nucleus 5 °Ti. It may be pointed out that the shell 

model calculation involved an independent variation of both 

the effective charges (e , e^) to fit the two maxima of I El 

^ 52 

for the nucleus 5 °Ti. The effective charges for Cr and 
54 Fe were, however, determined by varying only the isoscalar 
charge, \ (e + e n ) , so as to fit the lower-momentum transfer 

region of IF! 2 while the isovector charge, \ (e p - e n ) , was 

. 50 m , 

kept fixed to that used in TI* 
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111*3.2 The Nuclei 46 ' 48 Ti, 5 °cr and 56 Fe 

' " “I » ■"■■■ 

Fig. III. 3 compares the calculated and the experi- 

mental |F f for the 0 + -*■ 2^ transitions in the nuclei 

46 , 48 m , 50 , 56 7 

Ti, Cr and Fe. Iwamoto et al, had carried out 

sne 11-model calculations for tne nuclei ^'^Ti and 8 ^Cr; 

we have also shown tneir results here for the sake of 

comparison. However, no calculations on form factors has 

so far been reported for the nucleus 88 Fe. 

As in the calculations reported in the previous 
section, the effective charges have been used such that the 
proton effective charge is always greater than the neutron 
effective charge by an amount e. Here we find that the form 
factors resulting from the PHFB wave functions are in strikingly 
good overall agreement with the experiments provided a 
reasonable variation of the effective charges - with e n ranging 
between 0.5 and 0.7 — is allowed. Thus, the values for 
proton and neutron effective charges that we have amployed 
are'(e , e ) = (1.7e, 0»7e) for 48 Ti, (l»6e, 0*6e) for 46 Ti 

Jr ^ 

and Cr, and (1.5e, 0„5e) for the nucleus 6 Fe. with these 
effective cnarges, the maximum discrepancy is seen to occur 
in the case of the nucleus Fe for q f f > 1*4 fra” wnere 

the present calculation underestimates the magnitude of 
the second peak by about 18 percent. - 

In the. nuclei 46 ' 48 ^ and ^Cr, the form factors resulting 
from the PHFB wave functions show significant improvement over- 
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factors for the o + - 2+ TRANSIT, ONS IN 
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form factors, yield B(E2) and Q values -which are consis- 

2 + 

tent with the currently available experimental information 
for all the nuclei considered here.- 

III. 4 Conclusions 

We have presented here the calculation of Coulomb form 
factors for the 0 + -*■ 2^ transition in some doubly, even Ti, 

Cr and Fe isotopes in terms of the projected HFB wavaf unctions. 
In keen contrast with the results obtained in some recent 
shell model calculations involving f + f 7/2 P 3/2 confi 9- 
ur at ions , it is seen that the use of the self -•cons is tent 
projected HFB wave functions resulting from a sligntly modified 
version of the Kuo— Brown interaction for the 2p-lf snail 
leads not only to a significant overall improvement in the 
calculated form factors, it e Iso obviates the necessity of 
using widely different effective charges for various nuclei. 
However, the application of Distorted Wave 3orn Approximation 
(DWBA) instead of PWBA employed in the present worK, is likely 
to tesul-t in further improvements in the magnitudes of the 
calculated form factors around- the second maxima. 

' ' In view of the intractability of large-scale shall 
model calculations in the 2p-lf shell nuclei, the method 
employed in the present work is expected to provide a reliable 
alternative for obtaining a reasonably- accurate microscopic 
description of the observed data on C2 form factors. 


.O 
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CHAPTER IV 


INELASTIC ELECTRON SCATTERING FORM FACTORS 
INVOLVING THE NON-YRAST 2 + TRVELS 
IN THE NUCLEI 4d Tj, ^cr AND 54/56 Fe 


IV . 1 Introduction 

Some recent inelastic electron scattering experiments"'’ 4 
have provided valuable data on the 0 + -> 2 1 transitions in the 

_ V 

40 50 54 4- » 

nuclei ( Ti, Cr, Fe) and the 0 -+ 2^ transition in the 

nucleus ~^Fe. Although the qualitative features of the 

observed form factors associated with the 0 + -► 2* 3 transitions 

are quite similar to those of the observed form factors for 

the 0 + h. 2^ transitions discussed in the preceding Chapter/ 

2 

the magnitudes of IF| for the former transitions are smaller 
by an order of magnitude than those for the latter. 

In this 'Chapter we discuss an extension of the formalism 
presented earlier with a view to provide an interpretation 

of the available form factors involving the non-yrast 2 + 

, . 48 . 50 54/56 

levels in the nuclei Ti, Cr and Fe. 

Recently Iwamoto et al . 5 nave also studied the C2 form 
factors involving the second excited 2 r states in some Ti 
and Cr isotopes in terms of the wave functions resulting from 
empiric a 1 effective interactions operating in the restricted 

valence space consisting of the £ 7/2 + ^7/2 ^3/2 

tions. In their work Iwamoto et al . 5 employed the empirical 
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effective interactions given by Oda et ai . as well as by 
7 

Yokoyama et al t . As discussed in the preceding chapter/ the 
results for the low momentum-transfer (q ££<1.5 ^ region 


of I F ! (0 + - 2p were consistent with experiments; the results 

were also quite stable towards the choice of different effective 
interactions. However, the choice of the f^y^ + ^7/2 ^3/2 
space proved too restrictive to permit an unambiguous inter- 
pretation of the available data for the high momentum-transfer 
part of the o"*" -* 2^ transition, as well as the 0 -*■ 2^ transition 

involving 0.5 fm _1 <q f < 2.0 fra” 1 ? the calculated results 
displayed large variation with respect to the choice of the 
effective interactions. Furthermore, none of the two 
affective interactions yielded acceptable simultaneous agreement 


for the 0 


, , , + 

0 .* 27 as well as the CT - 2 n transitions, although 


relatively large values of the neutron effective charges 
e_ = l.Oe ard e = 1.2e for the nuclei 48 Ti and 5 °Cr, respecti- 


n 

vely 


n 


were emp loyed . 


As we have already seen, the yrast wave functions 

projected from the HIB intrinsic states of prolate snapes 

resulting from the modified version of the KB effective 

interaction 8 ' 5 for the full 2p-lf space provide a good descri- 

ption of the observed form factors for the 0 2^ transitions 

in several 2p—l£ shell nuclei • Here we show that a fairly 

2 

satisfactory interpretation of the available IF I data for 

48 50 54 

the 0 + - 2 + 2 transitions in the nuclei ( Ti, Cr, Fe) 
and the 0 + ^ 2* transition in the nucleus 56 Fe can also be 
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attempted in terms of reasonable values of the proton and 

neutron effective charges by describing the yrare and 

yr states in terms of the ob late variational solutions 

. 11 

resulting from the modified KB interaction . 

In section IV. 2 we discuss the details of the calcu— 

lational framework. In section IV. 3 we present the comparison 

, + 

of the calculated 0 2_ form factors with the experiments. 

/ o 

Finally section IV. 4 contains some concluding remarks, 

IV . 2 calcu lational Framework 

The details concerning the calculation of the 
inelastic scattering form factors for the 0 + -*■ 2^ transition 
in terms of the states with J = 0,2 projected from the same 
intrinsic HFB state (with K=0) have already been discussed 
in section III. 2. 2. However the states with J=0,2 involved 
in the calculation of the 0 + - 2 1, 3 form factors necessarily 

belong to two different intrinsic states; whereas the J=0 
belongs to the yrast band, the states J = 2 ^ ( 2^) belong to 
the yrare ( yr ) band. Further, as discussed later, the 
yrast and the non-yrast bands belong to the prolate and oblate 
intrinsic states, respectively, in the nuclei considered in 
this Chapter. 

The eigenstates of J projected from the HFB state 
! $ k ( 6)> can be written as 

1 JXL 6) > = /d^CQ) r(Q)IS (6»d0 » 

8 n 


(IV.l) 
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where R(Q) and D^(Q) are the rotation operator and the 
rotation, matrix, respectively. Here 5 denotes the type 
(prolate or oblate) of quadrupole deformation of the intrinsic * 
state. 

The intrinsic state ^.,,(6 )> can be written as 
(Eq. 11.44) 

'*K (a)> = & <^(6> + v”(6) C-)l0> , (IV.2) 


where the creation operator bl can be expressed 


as 


b 


1 

im 


^( 6 ) 




b'~ = 
im 


^i(6, 


(-D 


j-m i 


j-m 


(IV. 3) 


As before, here j labels the spherical single-particle orbits 
(lf 7/-2' 2p 3/2 / 2p l/2' lf 5/2* * 

. ,2 

The PWBA expression for the squared form factor IF I 

involved in the Coulomb scattering is In this case given 

. 12 
by 

lF ^> |2 = *5 nrar l<J f (6 2 )llf L ,|J i (6 i )>|2 > (IV - 4 > 

Z i 

wnere Z is the atomic number of the target nucleus. Here 
J i ( 6^ and J f ( <$ 2 ) denote the initial and the final scattering 
states resulting from angular momentum projection on the 
intrinsic states 1$ K ( 6^ > and I^ K (6 2 )>, respectively. The 
one-body operator is given by the eqn. ( III. 23) . 
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/s. 

The expression for the matrix element < J £ (6 ) If |J.(6 1 )> 

£ Z LM X -L 

is as follows: 


< J f (S 2> r V J i ( *l )> = t nJ (6i)n J '^2>]" 1/2 T? 

^ O f* 


J 2 J' 
-V V 0 


<o‘®> 


x n( 6 1 / 6 2 ,©)b 2 [ 2 < a|f p>P a p( 5 1 , 5 2 /®) ] sln©d© 

3 3 


(IV. 5) 


where 

n J (6) = J [det (1+F(6,9) f *(6 ) ) ] 1//2 d^(©)sin©d© , (IV. 6) 
n(6 l' 6 2 /0) = [det ( l+F(6 1/ ©)f T (6 2 )) ] 1/2 (IV. 7) 


and 

t t 

P ap ( W 01 =[M(6 1 ,6 2 ,e)/(l + M(6 1 ,6 2/ 9))] a 3 ^ 

(IV. 8 ) 

with 1*1(6^ 2 ,©> « F(6 1/ ©)f r (6 2 ) . (iv. 9 ) 


The matrices f 1 and F are defined by equations analogous 
to the ones given earlier (see Eqs.II.56 and 11.63): 


m 


m 


(6) - | c/ i( 6) 

CL p 


V. 


a 


yiti „ m /-m, 


(IV. 10) 


’i a a 


and 




.2 , d J<X ,(©) d~^ ,(©) f _(6). (IV. 11) 

m a m p m aPa m p m p 


ap 


The calculations were performed as follows. First the 
wavef unctions $ (6^ and $ Q (6 2 ) were obtained, for the yrast and 
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the yrare (or yr ) bands respectively, by carrying our the 

HFB (or the HE') calculations. The amplitudes (U^(6), V™(6)) 

and the expansion coefficients 0^.(6) were used to set up 

J **- " 

the £(6) matrices. Then F(<5,0), M( 6 , 6 ©) and 1/[1+M( 6^6 ,©) ] 

matrices were evaluated for twenty values of the Gaussian 
quadrature points between the range ( 0 , 11 / 2 ) . The form factor 
was then computed’ from equation (IV. 5). 

The centre-of-mass correction and the proton— finite- 
size effect have been taken into account by multiplying 
factors identical to the ones employed in the preceding 
Chapter. 

IV. 3 Results and Discussion 

Table TV. 1 summarizes our results concerning the 

variational solutions associated witn the (yrast, yrare) 

4t> 50 54 

levels in the nuclei ( Ti, cr, Fe) and the (yrast, yr ) 

5 6 

levels in the nucleus Fe. 

4b 5 4 

In the case of the nuclei Ti and Fe, the excited, 
self-consistent solutions of either the HFB or the HF-3CS 
type do not exist; the inclusion of pairing correlations 
always leads one to the variational solution belonging to 
the yrast band, irrespective of the choice of the starting guess 
for the wavef unctions. The non-yrast states have, therefore, 
been projected from the (excited) HF solutions in these nuclei. 
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Further, as mentioned in the preceding Chapter, a description 

54 

of the yrast states in Fe in terms of the HF-8CS (rather 
than HFB) solution is necessitated because of the shell 
closure at d= 28 . 

The suitability of the variational solutions considered 
here vis-k-vis the description of various yrast as well as 
non-yrast 2 + levels can be judged from the following consi- 
derations. It is reasonable to assume that the difference in 
the variational energies of two solutions represents roughly 
the difference in the energies of the 2 + states projected from 
these solutions, provided the two solutions have nearly the 
same absolute value of the quadrupole moment. From the 

results presented in column 6 and 7, Table IV. 1, we note 

48 

that the variational solutions obtained in the nuclei. Ti 

54 2 

and Fe satisfy this requirement - the 1<Q 0 > I values for 

the prolate and oblate solutions do not differ much. Thus 

the relative separations of the non-yrast and yrast 2 + states 

suggested by the intrinsic energies are about 1.4 MeV and 

54 

1.3 MeV, respectively, for the nuclei Ti and Fe. It 
is satisfying to observe that these estimates are quite 
close to the observed separations between the various 2 + 
states in these nuclei. 

The difference in the energy of an intrinsic state 
and that cf the 2 + state projected from It depends on the 
magnitude of deformation of the intrinsic state. In the 
case of the nuclei ~^Cr and ^Fe the intrinsic deformation of 
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the minimum-energy solutions are significantly larger than 
those of the excited variational solutions. This suggests a 
larger value of the difference [E H p B (Prolate)-E( J= 2 ^) ] 
compared to the difference [E^gC Oblate) -E( J =2 )] . 
Consequently, the separation between the yrare/yr 2 + state 
and the yrast 2 + state is expected to be somewhat larger 
than the difference in the corresponding intrinsic energies. 
Keeping this in mind, it is seen that the differences in 
the calculated variational energies for the nuclei 5 °Cr and 
54 

Fe are remarkably consistent with the observed relative 
energies of the 2 + states. 

In the form factor calculations presented here we 
have errployed the A- independent effective charges 
(e / e ) = (1.4e, 0.4e) for the transitions 0 + - it These 
values of effective charges are somewhat smaller than the 
ones errployed in our calculations of the 0 + -* 2 £ form factors 

presented in Chapter III. The reduction in the values of 
the effective charges in going from the yrast 2 + states to 
non-yrast 2 + states is not entirely unexpected since the 
states 2 ^ are usually more susceptible to core-polarization 
effects compared to the states it, 3 * 

2 4 

Figure IV. 1 compares the calculated and experimental ' 

IF 1 2 for the 0 + -> 2^ as well as 0 + -+ 2* transitions in the 
48 

nucleus Ti. We have also given here the results obtained 
by Iwamoto at al . 5 in the framework of the shell model involving 
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the £rj ^ 2 + ^ 7^2 ^3/2 confi 9 urations * In the earlier work. 

two sets of enpirical effective interactions were employed. 

One was the set of matrix element obtained by Oda et al «^ 

through a chi-square fit to the 38 levels chosen from Ca and 

Sc isotopes / and the other was the set obtained by Yokoyama 
7 

et al . through a fit to the 63 data selected from the N=27 

and 28 isotcnes (A = 47-55). Hereafter we shall refer to 

these sets of matrix elements as the v( Ca-Sc) and the 

V( 27-28) interactions, respectively. The broken and the dotted- 

dashed curves represent the results obtained by Iwamoto 

et al .^ with the enpirical effective interactions V( 27-28) 

and v( ca-Sc), respectively, in conjunction with the 
n n— 1 

f 7/2 + f 7/2 P 3/2 conf igurations. The solid curves give the 
results obtained in the present calculation. 

The results presented in Fig. IV. 1 bring out the sensi- 
tivity of the restricted shell model calculations for the nigh 
momentum- transfer part of the form factors associated with 
the 0 + -*■ transitions towards the choice of the effective 

1 r 2 

interactions. Further, none of the effective interactions 


employed earlier in conjunction with the f^^ + ^ 7/2 ^3/2 
space provides an adequate simultaneous interpretation of 
the available data involving both the 0 + -* 2* as well as the 

0 + 2^ transitions. Whereas the interaction V( Ca-Sc) leads to 

reasonable agreement with the experiments for q > 1.5 fm in 
the case of the 0 + •* 2* transition, the calculated values (not 
shown in the figure) for high moment um-transfer involving the 
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+ + 

0 -* 2 ^ transition are smaller by more than an order of 

2 

magnitude. The shell model estimates for I FI around its 
second maximum obtained with the interaction V( 27-28) , on 
the other hand, are smaller than the observed ones by roughly 
a factor of two for the 0 + -+ 2^ as well as the 0 + -»• 2"^ 
transitions. 

It is seen that the present calculation yields an 

adequate overall quantitative description of the available 

4 - + + + 

data for both the 0 -*• 2^ as well as the 0 -*■ 2 ^ transitions. 

In particular, the PHFB values for the squared form factor 
around q e ££ = 1.7 fm ^ show substantial irrprovements over the 
earlier results. Our choice of the neutron effective charge 
also contrasts keenly with the values employed in the earlier 
work; the present calculation involves the effective charge 
e n = °*^ e wherea s the restricted shell model calculation 
employed e n = l.Oe. The proton effective charge employed 
in the earlier calculation is, however, the sarre as the 
one used in the present one. 

A significant discrepancy between the calculated and 
the observed results occurs around the first maximum associ- 
ated with the 0 + -* 2* transition wnere the PHFB results 
overestimate the observed values by aoout a factor of two. 

In view of the sensitivity of the form factors for the 
0 + 2^ transition towards small configurational admixtures, 

an inclusion of the quasiparticle excitations in the oblate 
HF state is likely to result in further improvements. 
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We next present in Fig. IV. 2 a comparison of the 

5 50 

calculated and experimental results for the nucleus Cr. 

The shell model results obtained with various effective inter- 

48 

actions have also been presented* In the nucleus Ti only 

4 - + 

the part of the calculated form factor for the 0 -*■ 2 ^ 

transition Involving q > 1.5 fm ^ is found to be 

sensitive towards the choice of the effective interaction 

employed in the restricted shall model work* In contrast to 

this, even the low momentum- transfer part of the restricted 

shell model estimates for iFl^ (0 + -*- 2+) displays significant 

50 2 

interaction dependence in the nucleus Cr; the [Fj values 
resulting from V(ca-Sc) and V( 27-28) differ by more than a 
factor of 5 throughout the range 0.5 fm q eff < 2.0 fm \ 
Further, .both the interaction V( ca -Sc) as well as V( 27- 28) 
fail to reproduce even qualitatively the observed form 
factors for the 0 + -*• 2^ transition; the magnitude of iFl^ 
around its second minimum are much too low compared with the 
experiments. 

The use of the projected HFB wave functions for the 2* 
state leads to the required enhancement around the second 

2 

maximum. One also observes a shift of the dip in the If! 
values towards the lower moment urn- transfer region. Although 

— 1 

the PHFB results for the 0 + - 2* form factors for °* 8 

are still somewhat larger compared to the experiments, the 
present calculation is seen to yield significant, improvement 
vis— a— vis simultaneous description of the form factors for the 
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0+ ■* 2^ as well as the 0 + ■+ 2* transitions in the nucleus 3 °Cr, 

It may be mentioned here that the effective charges 

enployed in the earlier calculation, (e , e n ) = (1.&0/ 1.2s), 

were also considerably larger than the ones involved in the 

present calculation. In view of the relatively -.large values 

of the calculated intrinsic quadrupole moment (see column 4, 

Table IV. 1) , one expects an enhanced involvement of the 

configurations outside the ( f^^ + ^ 7/2 ^3/2^ s P ace ^ t *’ le 
, 50 

nucleus Cr, The increased value of the effective charges 

5 

employed in the earlier work may be intended to simulate 
this effect to some extent. 


Figure TV. 3 presents a comparison of the calculated 

and the experimental 3 ]Fj 2 for the 0 + - 2+ and 0 + -* 2^ transitions 
54 

in the nucleus Fe. Although the present calculation yields 
an adequate description of the form factors around the first 
maximum for both the 0 + -+ 2* as well as 0 + ■+ 2* transitions,, 
one notices significant discrepancies for JF{ 3 (0 + -» 2^ 
for q e£f > 1.4 fm 1 . In particular, the value of the form 
factor for 0 -*• 2 ^ transition around its second maximum is 

underestimated by roughly a factor of 6. In view of the 
smallness of the contribution of the neutrons towards the 
quadrupole moment of the oblate HF intrinsic state associated 
with the 2+ state (see column 4, Table IV. 1), the observed 
discrepancy may be due to insufficient admixture of various 
shell model. configurations involving neutrons. 



0.5 


1.0 1.5 2.0 

q eff (frn 1 ) 

FIG IV. 3 THE EXPERIMENTAL AND CALCULATED 
SQUARED FORM FACTORS FOR THE 
, TRANSITIONS IN 54 F C . 
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3 4 

Figure TV. 4 coirpares the calculated and experimental * 

form factors for the 0 + -* 2^ and 0 + -* 2* excitations in the 

, 56 

nucleus Fe. It turns out that the present calculation 
provides acceptable overall quantitative agreement with the 
available data simultaneously for the 0 -* 2^ as well as 

0 + •* 2* transitions. 

The experimental IfI^ data on 0 + -*■ 2^ transition for 
q > 1.4 fm” 1 is characterized by large error bars. Avail- 
ability of more accurate data for the momentum transfer range 
1.5 <q< 2.0 fm 1 will permit a better assessment of the 
present calculational framework;. 

IV. 4 Conclusions 

Summarizing, we have shown that the use of the yrast, 
yrare as well as yr J = 2 states projected from self- 
consistent prolate and oblate intrinsic states resulting from 
the modified KB interaction for the 2p-lf shell permits a 
reasonably adequate simultaneous description of the observed 
data involving the 0 + -* 2* 0 transitions in the nuclei 

48 Ti, 5 °Cr and 56 Fe which possess sizable quadrupole defor- 
mation. 

As mentioned earlier, an important feature that 
characterizes the observed form factors for the 0 + 2"^ 3 

transitions in the nuclei Ti, Cr and Fe is their 
reduced magnitude — by an order of magnitude compared to 
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those for the 0 + -* 2^ transitions. Present microscopic 
description offers a simple qualitative understanding of 
this feature of the 0 + -* 2"t _ transitions in terms of the 
reduced overlap between the initial (0 + ) and the final (2 
states due to their different intrinsic parentage. 


+ r 1 



REFERENCES 


1. J.W. Lightbody, Jr., et al ., Bull. Am. Phys. Soc. 20 , 

568 ( 1975) . 

2. A. Hotta at al ., Res. Rep. Nucl. Sc., Tohoku Univ. .9,7 

( 1976)7 10, 18 ( 1977). 

3. K. Hosoyama et al .. Res. Rep, Nucl. Sc., Tohoku Univ* 

11, 1 (1978). 

4. J. Heisenberg, J.S. McCarthy, and I. Sick, Nucl. Phys. 
A164, 353 ( 1971) . 

5. T. Iwamoto, H. Horie, and A. Yokoyama,. Phys. Rev. C25 , 

658 ( 1982) . 

6. T. Oda, K. Muto and H. Horie, Lett. Nuovo Cimento, 18, 

549 ( 1977) . 

7. A. Yokoyama, T. Oda and H. Horie, Prog. Theor. Phys. 60, 
427 ( 1978). 

8. T.T.S. Kuo and G.E. Brown, Nucl. Pnys. A114 , 241 ( 1968). 

9. J.B. McGrory, B.H. Wildenthal and E.C. Halbert, Phys. 

Rev. C2, 186 ( 1970). 

10. G. Mukherjee and S.K. Sharma, Phys. Rev. C29, 2101 ( 1984). 

11. G. Mukherjee and S.K. Sharma, Phys. Rev. C31 , 689 ( 1985). 

12. H. TJberall, Electron Scattering from Corrplex Nuclei 
(Academic, New York, 1971). 


4 



CHAPTER V 


CONCLUSIONS 


Inelastic electron scattering easily surpasses many 
other methods employed for investigating nuclear structure 
as far as versatility, possibility of unambiguous interpre- 
tation and accuracy is concerned. The nuclear matrix elements 
extracted from the measurement of the cross section for the 
inelastic Coulomb scattering can be considered to be the 
analogs of the usual matrix elements of the electric 2 —pole 
operator at non-zero momentum- transfers* In contrast to 
the ordinary Y -transitions, the momentum-transfer depen- 
dence of the nuclear matrix elements contains information 
about the spatial structure of the nuclear states involved 
and can therefore provide us a good test of the nuclear wave 
functions* 

The work presented in this thesis was motivated 
primarily by a desire to assess the validity of the 
projected HFB method in the context of a microscopic discu- 
ssion of the recent Coulomb form factor data involving the 
low-lying 2 + levels in some 2p-lf shell nuclei? such an 
assessment is warranted because one is constrained to use 
the projected HFB method as an alternative to large-scale 
shell model calculations which are intractable In most of 
the %>-lf shell nuclei. Our results indicate that the 



85 


ca lculat ional framework employed here can be considered 
reliable for interpreting and correlating the 
available data on the usual electromagnetic properties — 
e.g* the static moments and the E2 transition probabilities — 
as well as the elect roexcitation form factors involving 
low-lying 2 levels. 

The results presented in Chapter III have revealed 

noticeable discrepancies between the calculated and the 

2 —1 
observed (F(q) I values around the second maxima (q^l.7 fm ) 

associated with the 0 + -» 2^ transitions in some nuclei. In 

this connection it will be worthwhile to examine the 

application of distorted wave Bom approximation (DW3A) 

instead of the PwBA. 

The results presented in Chapter IV indicate that the 
projected HFB prescription leads to an overestimation of 
the 0 + -♦ 2* form factors around q^O.? fm - '*' in some cases 
although the effective charges involved in these transitions 
are slightly smaller than the ones used for the 0 + -=► 2^ 
transitions. This appears to necessitate an inproved 
description of the yrare 2 + levels via an inclusion of some 
additional quasiparticle-excited configurations. 



APPENDIX A 


THE ( lf 7/ / 2 A 2P3/2' 2pi/ 2 ' lf 5 / 2 ^ MATRIX ELEMENTS 
OF THE KB AS WELL AS THE' MWH INTERACTIONS 

We tabulate here the matrix elements <abJT I V I cd JT > 
of the Kuo-Brown (KB) interaction. The columns labelled 
by (ayb/C/d) contain twice the j-value of the various 
shell model orbits. 


T 

a 

b 

c 

d 

J 

< I V l> 

J 

< I V 1 > 

J 

< Iv !> 

1 

7 

7 

7 

7 

0 

-1.807 

2 

-0.785 

4 

-0.087 






6 

0.226 





1 

7 

7 

7 

5 

2 

0,000 

4 

-0.406 

6 

-0.716 

1 

7 

7 

7 

3 

2 

-0.502 

4 

-0.307 



1 

7 

7 

7 

1 

4 

-0.293 





1 

7 

7 

5 

5 

0 

-2.788 

2 

-0.638 

4 

-0.400 

1 

7 

7 

5 

3 

2 

0.305 

4 

0.305 



1 

7 

7 

5 

1 

2 

-0.470 





1 

7 

7 

3 

3 

0 

-0.783 

2 

-0.269 



1 

7 

7 

3 

1 

2 

-0, 250 





1 

7 

7 

1 

1 

0 

-0.714 





1 

7 

5 

7 

5 

1 

-0.287 

2 

-0.096 

3 

0.022 






4 

0.031 

5 

0.156 

6 

-0.894 

'l 

7 

5 

7 

3 

2 

0.014 

3 

-0.118 

4 

-0.114 






5 

-0.010 





1 

7 

5 

7 

1 

a 

-0.012 

4 

-0.093 
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T 

a 

b 

c 

d 

J 

< Iv | > 

J 

< IV 1> 

J 

<m> 

1 

7 

5 

5 

5 

2 

-0.631 

4 

-0.468 



1 

7 

5 

5 

3 

1 

-0.03S 

2 

0. 295 

3 

-0.132 






4 

0.493 





1 

7 

5 

5 

1 

2 

-0.342 

3 

0.080 



1 

7 

5 

3 

3 

2 

-0.018 





1 

7 

5 

3 

1 

1 

-0.079 

2 

-0.150 



1 

7 

3 

7 

3 

2 

-0.861 

3 

-0.027 

4 

-0.047 






5 

0.145 





1 

7 

3 

7 

1 

3 

0.078 

4 

-0.502 



1 

7 

3 

5 

5 

2 

-0.446 

4 

-0.156 



1 

7 

3 

5 

3 

2 

0.402 

3 

0.026 

4 

0.619 

1 

7 

3 

5 

1 

2 

-0.966 

3 

-0.059 



1 

7 

3 

3 

3 

2 

-0.325 





1 

7 

3 

3 

1 

2 

-0.319 





1 

7 

1 

7 

1 

3 

0.029 

4 

-0.274 



1 

7 

1 

5 

5 

4 

-0. 252 





1 

7 

1 

5 

3 

3 

0.140 

4 

0.531 



1 

7 

1 

5 

1 

3 

-0.078 





1 

5 

5 

5 

5 

0 

-0.860 

2 

-0.218 

4 

0.298 

1 

5 

5 

5 

3 

2 

0.023 

4 

0.102 



1 

5 

5 

5 

1 

2 

-0.186 





1 

5 

5 

3 

3 

0 

-0.777 

2 

-0.128 



1 

5 

5 

3 

1 

2 

-0. 256 





1 

5 

5 

1 

1 

0 

-0.392 
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T 

a 

b 

c 

d 

J 

<1Vi> 

J 

< lvl> 

J 

<!vl> 

1 

5 

3 

5 

3 

1 

-0.032 

2 

0.215 

3 

■ 0. 120 






4 

-0.243 





1 

5 

3 

5 

1 

2 

0.389 

3 

0*007 



1 

5 

3 

3 

3 

2 

0.068 





1 

5 

3 

3 

1 

1 

-0.058 

2 

0.159 



1 

5 

1 

5 

1 

2 

-0.135 

3 

0.205 



1 

5 

1 

3 

3 

2 

-0.047 





1 

5 

1 

3 

1 

2 

-0.241 





1 

3 

3 

3 

3 

0 

-1.206 

2 

-0.380 



1 

3 

3 

3 

1 

2 

-0.601 





1 

3 

3 

1 

1 

0 

-1.465 





1 

3 

1 

3 

1 

1 

0.152 

2 

-0.688 



1 

1 

1 

1 

1 

0 

-0*249 





0 

7 

7 

7 

7 

1 

-0.5 25 

3 

-0.208 

5 

-0.502 






7 

-2.199 





0 

7 

7 

7 

5 

1 

1.894 

3 

1.005 

5 

0.-901 

0 

7 

7 

7 

3 

3 

-0.482 

5 

-0.816 



0 

7 

7 

7 

1 

3 

0.640 





0 

7 

7 

5 

5 

1 

1.071 

3 

0.517 

5 

0.170 

0 

7 

7 

5 

3 

1 

0. 227 

3 

-0.087 



0 

7 

7 

5 

1 

3 

0.033 





0 

7 

7 

3 

3 

1 

-0. 279 

3 

-0.296 



0. 

7 

7 

3 

1 

1 

0.392 





0 

7 

7 

1 

1 

1 

0.184 
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T 

a 

b 

c 

d 

J 

< Iv j> 

J 

< I V 1 > 

J 

<fv|> 

0 

7 

• 

5 

7 

5 

1 

-3.621 

2 

-2.731 

3 

-0.985 






4 

-1.886 

5 

-0.112 

6 

-2.217 

0 

7 

5 

7 

3 

2 

-0.735 

3 

0.3 28 

4 

-0.039 






5 

0* 405 





0 

7 

5 

7 

1 

3. 

-0-. 440 

4 

-0.674 



0 

7 

5 

5 

5 

1 # 

-0.416 

3 

0.675 

5 

1.138 

0 

7 

5 

5 

3 

1 

0.953 

2 

-0.705 

3 

0.628 






4 ’ 

-0.544 





0 

7 

5 

5 

1 

2 ‘ 

0.427- 

J ‘ ' 

0.539 



0 

7 

5 

3 

3 

"1 ' 

0.871 

3 

0.538 



0 

7 

5 

3 

1 

’I ' 

-1.449 

2 ‘ 

-0.659 



0 

7 

5 

1 

1 

1 

- 0*168 





0 

7 

3 

7 

3 

2 * 

-0.'293 

3 

-0.604 

4 

-0.164 






5 ‘ 

-2.165 

> 




0 

7 

3 

7 

1 

3 _ 

1. 295 

4 

0.108 



0 

7 

3 

5 

5 

3 _ 

0.158 

5 

0.048 



0 

7 

3 

5 

3 

2_ 

-1.123 

3 

-0 . 408 

4 

-0.614 

0 

7 

3 

5 

1 

2 

0.838 

3 

0.518 



0 

7 

3 

3 

3 

3 ' 

-0.511 

• 


... • ■ 


0 

7 

3 

3 

1 

l' 

-0.481 





0 

7 

1 

7 

1 

3 ' 

-1.484 

4' 

-0.746 



0 

7 

1 

5 

5 

3 \ 

—0. 270 

* 















0 

7' 

1 

5 

3 

3 ' 

0.139 

* 

4 

-1.188 



0 

7 

1 

5 

1 

n 

3, 

0.400 



* 


0 

7 

1 

3 

3 

3 

0.461 
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T 

a 

b 

c 

d 

J 

< lvl> 

J 

<Iv }> 

J 

<lv 1> 












0 

5 

5 

5 

5 

1 

—0*01 2 

3 

-0.196 

5 

-1.687 

0 

5 

5 

5 

3 

1 

-0. 29 2 

3 

-0.327 



0 

5 

5 

5 

1 

3 

-0.391 





0 

5 

5 

3 

3 

1 

0.040 

3 

-0.182 



0 

5 

5 

3 

1 

1 

0. 159 





0 

5 

5 

1 

1 

1 

-0.084 





0 

5 

3 

5 

3 

1 

-2.164 

2 

-1.236 

3 

-0.459 






4 

-0.722 





0 

5 

3 

5 

1 

2 

0.255 

3 

-0.988 



0 

5 

3 

3 

3 

1 

-0.073 

3 

-0.340 



0 

5 

3 

3 

1 

1 

0.620 

2 

-0.071 



0 

5 

3 

1 

1 

1 

-0.393 





0 

5 

1 

5 

1 

2 

-0.081 

3 

-1.234 



0 

5 

1 

3 

3 

3 

-0.084 





0 

5 

1 

3 

1 

2 

0.464 





0 

3 

3 

3 

3 

1 

-0.639 

3 

-1.832 



0 

3 

3 

3 

1 

1 

1.554 





0 

3 

3 

1 

1 

1 

0.709 





0 

3 

1 

3 

1 

1 

-2.291 

2 

-1.897 



0 

3 

1 

1 

1 

1 

0.686 





0 

1 

1 

1 

1 

1 

-1.074 
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The MVfH interaction referred to in the thesis is 
obtained from the KB interaction by 1 substituting the 
following matrix elements in place of their counterparts. 


T 

a 

b 

c 

d 

J 

< fv i > 

J 

< I v| > 

J 

< 1 V 1 > 

1 

7 

7 

7 

7 

0 

-2.107 

2 

-1.115 

4 

-0,097 

1 

7 

5 

7 

5 

_L 

-0.037 

2 

0.154 

3 

0,27 2 






4 

0.281 

5 

0.406 

6 

-0,644 

1 

7 

3 

7 

3 

2 

-0.561 

3 

0.253 

4 

0,283 






5 

0.485 





1 

7 

1 

7 

1 

3 

0.279 

4 

-0.0 24 





